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Abstract 

We study rj c and J/tp hadroproduction induced by multiple scattering 
off fixed centres in the target. We determine the minimum number of 
hard scatterings required and show that additional soft scatterings may 
be factorized, at the level of the production amplitude for the r] c and of 
the cross section for the J/ip . The «//?/> provides an interesting example 
of soft rescattering effects occurring inside a hard vertex. We also explain 
the qualitative difference between the transverse momentum broadening 
of the J/%b and of the T observed in collisions on nuclei. We point out 
that rescattering from spectators produced by beam and target parton 
evolution may have important effects in J/ij; production. 


1 Work supported in part by the EU/TMR contract ERB FMRX-CT96-0008. 



1 Introduction 


Quarkonium production is a sensitive measure of soft rescattering effects in hard 
collisions. Whereas the creation of a heavy QQ pair is a process of scale ttiq, 
the binding energy of the heavy quarks in non-relativistic quarkonium is only 
on the order of a 2 mg. Thus for charmonium the binding energy is of order 
2m d — wjm — 600 MeV, which is only moderately larger than A qcd- 


There is considerable debate concerning the correct theoretical description of 
quarkonium production. Whereas the ‘color singlet mechanism’ (CSM) [pj ap¬ 
pears to give good agreement with J/i[j photoproduction data 0, it fails by more 
than one order of magnitude in hadroproduction || 0. |5|. |J. This has motivated 
proposals of production mechanisms where radiation at the binding energy scale 
plays an essential role, such as ‘color evaporation ’0BB0 and the ‘color octet 
mechanism’ [111] , [T^ . p~3|] . All models seem to have problems with some aspects of 
the data, however. 


The observed nuclear target A-dependence |T4| of the quarkonium cross section 
gives further insight into the production process. Compared to lepton pair pro¬ 
duction at the same hardness, quarkonium production shows much bigger nuclear 
effects. The cross sections (averaged over xf ~ 0) can be parametrized as 


a A - cr N A c 


with 


a ~ 1.00 
a ~ 0.92 ± 0.008 
a ~ 0.962 ±0.014 


p/p 

T(IS') 


( 1 . 1 ) 


These values are taken from Refs. ||15|, pR , |T/j |. respectively. 

Furthermore, there is evidence that the p_i_-broadening induced by the nucleus 
depends on the quark mass [IT]. [TS| , 

0.113 ±0.016 GeV 2 


(pi (A)) - (pi( z H)) = { 0.34 ± 0.08 GeV 2 

0.667 ±0.133 GeV 2 


p 1 p 

T(IS') 


( 1 . 2 ) 


where A = 184 (W) for the E772 data on p + p and T production, whereas 
A = 195 (Pt) for the NA3 data on J/if production. 


The sizeable nuclear effects in quarkonium production shown by Eqs. (|1 . 1|) and 
(|1 . 2|) raise two questions, which we shall consider in this paper. 


(i) Can 3 Si quarkonia be produced through a hard gg —> QQ subprocess, with 
an additional soft rescattering gluon coupled to the QQ pair? If not, how 
does soft rescattering factorize from the hard production amplitude? 

(ii) A color octet heavy QQ pair of high momentum should behave like a point¬ 
like gluon in soft rescattering processes. How can the p^-broadening depend 
on the quark mass, as indicated by Eq. o? 
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These questions concern the interplay of soft rescattering and hard production 
processes, which have so far not been thoroughly studied in QCD |0|. Soft 
rescattering by itself has been studied in the problems of high energy parton 
propagation and energy loss in a dense or hot medium | 2 (| . 


Hard production amplitudes have been calculated at leading twist, involving only 
a single scatterer in the target. A factorization between hard and soft processes 
a is usually assumed to hold also in quarkonium production - which is indeed 


necessary for any firm perturbative predictions. The puzzles of quarkonium pro¬ 
duction and the experimental evidence for rescattering effects that depend on the 
hard scale motivate a more detailed investigation. 


We shall study the rescattering and production of quarkonium in the limit of 
(asymptotically) high energy. Since the hardness of the process (as measured 
by the quark mass) is kept fixed, this implies that the momentum transfer to 
the target is a (vanishingly) small fraction of the projectile energy. Such a limit 
seems natural for much of the data, which is at high energy and appears to obey 
Feynman scaling. 


The high energy limit considerably simplifies elementary cross sections. For ex¬ 
ample, it is easy to verify that the e/i —> e/z cross section reduces to the Ruther¬ 
ford one when the electron momentum tends to infinity while the transverse 
momentum transfer and the muon momentum are kept fixed. The fact that the 
scattering cross section in this limit is independent of the target momentum and 
mass appears to be quite general, and holds at least to lowest order for elementary 
targets. 


We shall take advantage of this simplification by modelling the target partons 
by very heavy quarks. This suppresses energy transfer and selects Coulomb 
exchange, thus considerably simplifying the calculation. In the Appendix we 
show explicitly (using the 7 g —> qq process) how the number of heavy quarks 
in our kinematic limit is related to the gluon distribution G(x —> 0). A similar 
approach has earlier been used for describing deep inelastic lepton scattering in 
the target rest frame and the creation of rapidity gaps [| 22 | . 

Rescattering effects in quarkonium and Drell-Yan production have been consi¬ 
dered previously, see, eg , Ref. |23|. Our assumption of a simple target structure 


allows a systematic and precise investigation of the multiple scattering amplitudes 
in QCD. This reveals interesting aspects of factorization between hard and soft 
processes, color dynamics and /^-broadening. 


We consider processes of the type shown in Fig. [I[. An incoming gluon of asymp¬ 
totically high energy mutiply scatters in the target producing a heavy quarko¬ 
nium. We use charm to represent a typical heavy quark and consider both 1 So 
(r] c ) and 3 Si (J/i/j) production. We assume that the bound state is produced in 
a color singlet state, ie, the cc pair couples directly to the charmonium through 
its wave function at the origin. Thus we do not consider higher order terms in an 
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expansion in the relative velocity of the charm quarks, which are relevant for the 
color octet mechanism 
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Since we do not consider gluon radiation, at 
least two target scatterings are required by charge conjugation invariance to pro¬ 
duce a while one is sufficient for the rj c . We are interested in understanding 
the systematics between hard and soft Coulomb gluon exchanges. 


q Pi 

-*- *—>■ 

i-<-i- - -«- 

q ■ ; -Pi 
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Figure 1: General amplitude for q c or J/ijj production induced by n scatterings 
off static centres 1,... n located at xi,... x n . The static centres are ordered and 
labelled according to their increasing longitudinal position, and ki is by definition 
the momentum transfered to the qq pair by centre # The color indices of the 
incident gluon and of the exchanged gluons are denoted by c and a;, respectively. 

A QED illustration will define more precisely what we mean by ‘hard’ and ‘soft’ 
momentum transfers in this work. The Rutherford scattering e + A —> e + A of 
a charged particle off a heavy nucleus is described by the Feynman diagram of 
Fig. |^a. Due to the Coulomb photon propagator the amplitude is proportional to 
1 /kj_, where k± is the transverse momentum transfer. The cross section is then 
given by an infrared divergent integral, a oc / d 2 k±/k^. For scattering on neutral 
atoms, the infrared cutoff is given by the inverse atomic radius R, 

k± is soft k± ~ — (1.3) 

R 

On the other hand, in the pair creation process of Fig. |2:b, 7 + A —» e + e _ + A, the 
scattering occurs off the neutral e + e~ pair. The amplitude is now proportional 
to the dipole moment oc l/m e of the pair, which introduces a factor k±/m e . The 
integral over momentum transfers is logarithmic, cr oc. J d 2 k±/k\, and has support 
from an extended region, 


C I 1 


11c, J/V 


p - 2 P 1 


k± is hard - « b « m e (1.4) 

it 

The electron mass sets the upper limit on the logarithmic integral. Indeed, for 
k± m e the virtual electron or positron in Fig. S> is far off-shell, and the 
fermion propagator causes the integral to converge. 
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Figure 2: (a) Rutherford scattering of a charged particle off a nucleus, (b) Pair 
creation process 7 + A — > e + e _ + A. 

The general classification of scattering into soft (monopole) and hard (dipole) 
exchange is not altered by effects due to a running of the coupling constant. In 
the following QCD analysis we shall thus for clarity assume a fixed coupling a s . 

In the QCD process of Fig. |lj. the incoming gluon has color charge, while the 
final charmonium is a color singlet. I 11 analogy to the QED case above we may 
expect that the gluon will multiply Rutherford scatter, with the infrared cutoff 
provided by the inverse nucleon radius Rr 1 ~ A qcd- At asymptotically large 
gluon energies the virtual cc pair travels a long distance, and it becomes highly 
probable that the pair is created far upstream of the target. Hence the multiple 
scattering in the target will actually occur off the quark pair rather than off the 
gluon. The compact pair is, however, created in a color octet state and it will 
behave like a gluon in soft rescattering processes. 

On the other hand, it is intuitively plausible that a scattering which changes the 
color state of the pair (from singlet to octet or vice versa ) will depend on the 
differing structure of quark pairs and gluons, ie, such a scattering will probe the 
dipole moment of the pair. As a consequence this scattering must be liardQ This 
applies, in particular, to the last scattering in Fig. ]Ij. which by definition changes 
the color state of the pair from octet to singlet. 

1 Conversely, a hard scattering may or may not change the qq color state. If the pair is 
turned into a singlet then any next scattering must also be hard, which implies a subleading 
contribution in the multiple scattering production process we are considering. 
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In this paper we shall show explicitly how the above expectations are realized in 
QCD, and also find the answer to the questions mentioned in the beginning. 

(i) Even though charge conjugation symmetry would allow the J/ip to be pro¬ 
duced through one hard and one soft scattering in the target (without 
gluon emission), this would break the factorization between hard and soft 
processes. The soft exchange cannot be reliably calculated and such a 
contribution would in fact make it impossible to predict J/ij) production in 
PQCD. We show that this problem does not arise since two hard scatterings 
are needed to produce a J/ip. 

(ii) We shall also see that the broadening of (p 2 ± ) in nuclei does in fact depend 
on the hard scale, as suggested by the data of Eq. Q- Multiplying the 
differential cross section by a factor p\ makes the integral sensitive to the 
upper cut-off, and (p\) then depends (logarithmically) on the quark mass. 

2 r] c production 

2.1 Single gluon exchange 

The 1 So cc state p c has positive charge conjugation and thus couples to two gluons 
- the projectile and one Coulomb gluon from the target as shown in Fig. [|. 



Figure 3: Amplitude for p c production induced by one scattering. 
The amplitude for scattering from a fixed centre at position X\ is 


M{gg Vc) 


Rod 3 /' d 3 ki 


f ’ 1 (2vr) 3 <5 3 (g + h - 2pi) 


exp (—iki ■ xi) 




x (Mi + M2) 



Tr 7 0 


7°(^ 1 - $ 1 + m)tfV S : 


=0 


M 2 = Tr (T c T ai ) 


(pi — k \) 2 — m 2 + ie 


( 2 . 1 ) 
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Here the fixed centres are treated as heavy quarks of initial and final color charge 
Ai, A\. In the nonrelativistic limit of the charmonium bound state the momenta 
of the c and c are the same, p± = p-2, hence 


q + k = 2pi = P 


( 2 . 2 ) 


where k is the total momentum transfer to the target ( k = k\ in Fig. [5j). The 
charm quark mass is denoted by m and we have used the operator ||24|j 

Vs=0 = Tf 75 ^ 1 + m ^ ( 2 - 3 ) 

to project the cc pair onto the r/ c wave function at the origin, which is given by 
Rq. The incoming gluon polarization vector is denoted by e. 

The denominator in Eq. ( |2.1| ) is 

(pi — k\) 2 — m 2 = (q — pi) 2 — m 2 = — 2pi ■ q = —q ■ k (2.4) 


In the limit q° —> oo and allowing the incoming gluon to have a finite transverse 
momentum, ie, q ~ (q°, q±, q° — q\/2q°), we can write 


(pi — ki )' 2 — m 2 

k z 

m 2 . 


~ — 2m 2 L 

2m 2 , + q± ■ k± 

___ 

qO 

= m 2 + k \/4 


(2.5) 

( 2 . 6 ) 
(2.7) 


and obtain an expression for the iq c production amplitude whose kinematics de¬ 
pends only on k\ (since k — k\ here), 


M{g(q)g{k 1 ) -»■ rj c (2pi)) ^ iq° 


Ro g 3 


m- 


\J~2A 


Tim 


n~ic 

1 A’ A 1 


exp (—iki± ■ Xi±) 


k i± 


e 0[iv3 k l 


(2.S) 

where £^ pa is the fully antisymmetric Levi-Civita tensor. In the following, the 
direction of the incoming gluon will be chosen along the £-axis, ie, we take qj_ = 0 
from now on. 


We can now make the following observations. 


• The rj c production amplitude A i(gg —> r/ c ) is proportional to ki±/kf ± , ie, 
the effective values of k\_ are hard, 

Rqcd “C ki± m (2-9) 

which justifies the use of perturbative QCD for this process. 
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• The fact that the effective upper limit of k\± is given by m originates from 
the quark propagator, proportional to l/m\ = l/(m 2 + fc 2 _|_/4). As we shall 
see and as already announced in Eqs. (|1.3| ) and (|1.4|) , the dependence of 
m _l may be safely neglected when calculating amplitudes or cross sections. 
However, this dependence can be important in the calculation of other 
quantities, like (p 2 ^) for g c or J/Vq since weighting the differential cross 
section by p\ can shift the effective values of the transverse momenta to 
be of order m. For this reason we will keep the m± dependence in the 
following. 

• In this process, the incoming cc pair is in a color octet state, but ki± is not 
soft as it would have been in the case of elastic gluon scattering. In order 
to turn the color octet cc into a singlet the exchanged gluon has to probe 
the color dipole moment of the pair. 

In the next subsection we consider the case when the cc pair undergoes two 
Coulomb scatterings in the target. 

2.2 T] c production through two gluon exchange 

According to the notation of Fig. [l] we assume momentum transfers k \, k 2 
from the two scattering centres located at x\ , x 2 , respectively. The longitudinal 
coordinates of the centres are taken to be x\ = 0, x\ = r > 0 (we do not consider 
double scattering from a single centre in this paper). The calculation, which is 
summarized below, shows the following. 

• ki± is soft, whereas k 2d is hard (in the sense of Eqs. (ED and Q) , ie, 

ki_i_ ~ A qcd] A qcd •C k 2 ± <C m ( 2 - 10 ) 

• The amplitude factorizes. It is given by a convolution of the incoming 
gluon elastic scattering amplitude gg —► g and the g c production amplitude 
99 -»• Vc given by Eq. ( pf ), 

r d 2 k\± 1 

iM{ggg -> g c ) = J2 j j^y^eilg^ A )g{k 1 ) -»• g(q + fci, A')] ^ 

x iM{g(q +k!,\')g(k 2 ) -*• g c (q + k)} (2.11) 

where X 1 is the polarization of the intermediate gluon and 

iM e i[g(q, X)g(kf) -»• g(q + h, A')] = 2q°g 2 e(X) ■ e(X')f caid T^ M 

exp(— iki± ■ xi±) 




This equation may be diagrammatically represented as 




* 

i 



X' 




h k 2 

* 

2 


(2.13) 


Let us now present the derivation of Eq. ([2.11|) , which within our approach is 
exact in the limit q° —> oo. We write the amplitude as 


M(ggg^rj c ) = 


Rog 5 


rpa\ rpCL2 

zl A' A. 1 A' . 


X 


" d 3 k 1 

exp(-zfci_L • x x _ - ik 2L ■ x 2 ±) 
P P 


exp [irikl - k z )\ 

(A + B) (2.14) 


where the amplitudes A ( B ) correspond to the 4 ( 8 ) diagrams shown in Fig. [|a 
(b). 

Since we take r > 0 the integral over kf can be evaluated by closing the integration 
contour in the upper half plane. 

It is rather straightforward to see that none of the diagrams of Fig. [i]b contribute 
in our limit, ie, B = 0. Namely, 

- Since the Coulomb exchange is instantaneous, the diagrams B x ... B 4 that 
have crossed exchanges involve the creation of qq pairs with energies oc 
q° —> 00 . All these diagrams have quark propagators whose residues in the 
Irn k\ > 0 plane vanish asymptotically. 

- Diagrams B 5 and B e which involve the three-gluon vertex are also negligible. 
Onr static scattering centres and the q° —> 00 limit restrict the Lorentz 
indices to be g = v = 0 for the two gluons attached to the centres and 
cr = 0, 3 for the gluon attached to the qq. Hence the vertex is at most on 
the order of the asymptotically vanishing longitudinal momentum transfers 
(c/Eqs. (|2.6| ) and ( j2.23|) below). 

- In diagrams B 7 and B s the first rescattering gluon is attached to the pro¬ 
jectile gluon. In the q° —> 00 limit the g —> qq fluctuation will occur long 
before the target, and thus such scattering should be suppressed. Techni¬ 
cally this appears in the following way. The product of the gluon and quark 
propagators in diagrams By and B$ involves the denominator 

[(q + h ) 2 + ie] [(pi - k 2 ) 2 - m 2 + ie] ~ 

h± -k 2 ± 


2(<? 0 ) 2 [kl + ^-iem-k 


qv 


— le 
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and thus has two poles in the upper k{ half-plane, whose residues can easily 
be seen to cancel in the q° —> oo limit. 


Thus, only the four diagrams of Fig. |4ja are nonvanishing and contribute to the 
^4-term in Eq. (|2.14| ). Separating the color factors c tJ , the Lorentz traces Try,- and 


the denominators A^ we write 


where 


a = E^« = E 


hJ 


Tr 

ij A-■ 

i,j 


c n = C 12 = Tr(T c T ai T a2 ) 
c 22 = c 21 = Tr(T c T a 2 T ai ) 


(2.15) 


(2.16) 


Tr n = 

Tr 


+ ^ 

+ + m) 7°(-^ 1 + lf , 2 + m 

Tr 22 = 

Tr 

7°(|* 1 

~h 

+ m) 7 °(^! 

-ty 1 — ty 2 + m)rf 

$ 

to 

Tr 

7 °(^ 

~h 

+ m )i(~Pi 

+ ^ + 771)7° (^ 

Tr 2 i = 

Tr 

7 ° (ft 1 

~h 

+ m)f(-p i 

+ + m) 7°(^ 1 

A n 

= 

A 22 = 

[(Pi 

— /c ) 2 — m 2 

+ ie][(p 1 - k 2 ) 2 - 

A12 

= 

A 2 i = 

[(Pi 

— Aq ) 2 — nr' 

+ ie][(pi - k 2 ) 2 


m + ie ] 

Reversing the order of the 7 -matrices in Tr 22 and Tr 21 one easily finds 

Trn = — Tr 22 ; Tr 12 = -Tr 21 


(2.17) 


(2.18) 


(2.19) 


In QED Eq. fl2.19|) ensures a vanishing transition amplitude between a three- 
photon state (of negative charge conjugation, C = —) and a 1 S 0 positronium 
state (C = +). Due to the color charges of QCD the analogous amplitude does 
not vanish but is proportional to the fully antisymmetric structure constants f a bc 
of SU(3), 


j _^ _ [Tr(T c T“iT a2 ) Tr(T c T a2 T ai )] 


Tr 


11 


A 


+ 


Tr 


12 


11 


A 


12 j 


2 /caia2 


Tr n Tr 12 


.An A12 j 

An evaluation of the traces gives, in the q° —> 00 limit, 

Trn = — 2 i(q 0 ) 2 £0^3 £^ (ki + k 2 ) u 

Tri 2 = —2i(q 0 ) 2 £0^3 £^ (—k\ + k 2 ) 1 


( 2 . 20 ) 


( 2 . 21 ) 
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The denominators may be simplified using 


\2 2 , • _ „0/,z , ^1-L ' ^2_L 


(pi — k\) — m + ie ~ g u (k{ + 


+ ie) 


(Pi — k 2 ) 2 — m 2 + ie ~ —g° — k z — 


h± ■ k 2 _l 
g° 


— «e 


( 2 . 22 ) 


Integrating over k{ in Eq. ([2.14 ) by closing the contour in the upper half-plane 
picks according to Eq. ( |2.22| ) the pole 


K = k* + + j£ 

gU 


(2.23) 


such that 
dk{ 


/ Cl 

A exp[ir(kl - k z )} = i(q 0 ) 2 f caia 2 e 0 ^ 3 ^ 


x 


+ 


(—+ k 2 y 


[k\ + /C2 

q°k z (-q°) ' (- g 0)(g0^ + 2fc ix -fc 2 ±) 


(2.24) 


Note that the poles of the Coulomb gluon propagators lead to a subleading con¬ 
tribution in the g 0 —>► oo limit. Using Eq. (|2.6| ) for gl = 0 and 


\h± ■ k 22 _ | -C m]_ 

we find for the full amplitude of Eq. (|2.14[) 

Rod 5 


(2.25) 


M(ggg^g c ) = 


-iq 


X 


m 2 ± v2Airm 
d 2 ku 


rpai rjia2 £ - -H 

1 A^Ai A' 2 A 2 J ca i a 2 L 0fJ,u3 

1 


(2tt) 


exp(-z/Ci_L • Xix - i/C 2 _L ■ x 2 ±) 


r .2 p /2 

^ 1 ± ^ 2 ± 


Writing 


^(A) = -E[e(A)- e (V) 1 ^(V) 


(2.26) 


(2.27) 


A' 


directly leads to the factorized expression of Eq. (|2.11|) . 

Let us note that Eq. (|2.26| ) involves the denominator m\ = m 2 + [k\±_ + k 2 ±) 2 /4:, 
instead of (m 2 + k 2 ± /4) for the rj c production amplitude contained in Eq. ( |2.11| ). 
By calculating M.(ggg —> rj c ) as a function of k± and x 2 ± — xi±, one shows that 
the effective values of the transverse momenta satisfy k±± -C k 2 ± ~ k± for all 


\x 2 ± — xu_| (provided fcj_ A qcd), so that Eq. (|2.11|) follows. However, as 
mentioned earlier the k±± dependence of m± can be important when calculating 
(pj_), for which the correct expression to start with will be Eq. (|2.26p. 
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The k 2 ± factor in Eq. ( |2.26| ) explicitly shows that the second exchange is hard in 
the sense of Eq. (|1.4| ) . The amplitude to turn the color octet cc pair into a color 
singlet is thus proportional to the dipole moment of the pair. In the next section 
we study the analogous situation in J/-0 production, which differs from the case 
of the rj c since a minimum of two exchanges with the target is now required. 


3 J/'il' production 

3.1 Two gluon exchange 

The J/-0 is a 3 Si charmonium state with negative charge conjugation, and hence 
couples to a minimum of three gluons. The lowest order amplitude (which does 
not involve gluons in the final state) is thus given by the diagrams of Fig. |a, with 
the rj c replaced by the J/^. For the reasons discussed in the previous section the 
diagrams B\ and B 2 of Fig. [|b do not contribute to high energy J/ip production. 
Hence the amplitude may again be written as in Eq. ( |2.14|) with B = 0. In 
the expression for A, the only difference is that the projection operator (|2.3|) is 
replaced by a projection onto a vector state PJ}] , 

Ts=i = +m ) (3.1) 

where e(S z ) is the J/rp polarization vector, defined in the J/ip rest frame as 

e{S z — ±1) = ex = (0,1, ±z, 0)/\/2 for transverse J/^ 

e(S z = 0) = eL = (0,0j_,l) for longitudinal Jj\b (3.2) 

We shall use the Gottfried-Jackson frame, where the z-axis in the J/^j rest frame 
is taken parallel to the projectile momentum q. The invariant product e • w for 
any vector w can then be expressed, in the high energy limit we consider, as 


ex ■ w = — e± ■ w± 

q ■ w pi ■ w 

eL ■ w = m - 

q ■ pi m 


(3.3) 


From the trace expressions corresponding to Eq. (|2.17| ) it can be seen that for 
the J/ij; amplitude 

Tr n = Tr 22 ; Tr, 2 = Tr 21 (3.4) 


which differ from the p c case of Eq. ( |2.19| ) by a sign. It is in fact easily seen that 
Eq. (|3.4|) for the J/ij} generalizes in the case of n scatterings to 


Tr A = (-l)"Tr A * 


(3.5) 
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where A and A* are charge conjugated diagrams. Eq. ( |3.5| ) is readily understood 
as a consequence of charge conjugation invariance in QED. 

The expression for the amplitude A in Eq. (|2.14|) for J/ij} production is then 


A = [Tr(T c T ai T a2 ) + Tr(T c T a2 T ai )] 


Tr 


11 


A 


+ 


Tr 


12 


11 


A 


12 J 


- . dr 


Tth Tr 12 


M 2 J 


.A n Ai 

In the limit q° —> oo the explicit expressions for the traces are 

Tr n = -TV 12 = -4 m(q 0 ) 2 f3r >L (\) 


(3.6) 


(3.7) 


where (3 depends on the polarization (A) of the incoming gluon and on that of 
the J/tl) (T,L), 


/?t(A) — —e(A) • ex 

Pl(X) = 2e(\) ■ Pim/m 2 ± 


(3.8) 


For the general case of n scatterings (c/Fig. |l|) the trace for a given diagram is 

Tr = —4m(g°)7MA)(—l) n? (3-9) 

where rig is the number of Coulomb gluons attached to the antiquark. 

Using now Eqs. (|3.6| ) and (|3 . 71 ) in the expression (|2.14|) and performing the 
hf-integral by picking up the pole (|2.23|) we hnd (once again \ku_ ■ k 2 ±\ -C m\), 


M(g(X)gg -> J/ip) = -iq°m 

r d 2 k\± 

X J J2^ 


Ro g° 


ml V2Att7ti 


T a 4 A T% 2A2 d caia2 (3 T , L ( A) 


exp(-ik 1± ■ fu. - ik 2 ± ■ x 2 ±) ^ A 2± (3.10) 


Since the numerator of the J/rj} production amplitude (|3.10|) is proportional to 
feu. and k 2 ± it follows that both transferred momenta are hard in the sense of Eq. 
(|1.4| ). This contribution to J/ij) production is then a higher twist effect compared 
to the standard one involving a single target scattering and a radiated gluon, as 
is required by factorization between hard and soft processes in PQCD. 


Next we shall investigate the case of one ‘extra’ target scattering. This is of 
importance for understanding the general systematics of rescattering effects in 
J/ij} production and the p^-broadening effects in nuclei. 
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3.2 Three gluon exchange in J/'ijj production 

In the high energy limit the J/-0 production amplitude induced by three scatter¬ 
ings on static centres is given by the eight diagrams of Fig. [5|. 

From the p c result ([2.26|) we may expect that the last transfer (k 3 ), which turns 
the cc pair into a color singlet, will be hard. It is less obvious if either one of the 
first two exchanges can be soft. As already noticed, the propagating quark pair 
stays dominantly in a color octet state even when undergoing a hard scattering 
(except for the last one). Therefore, one would expect it to suffer soft rescattering 
both before and after the first hard exchange. This is indeed what we shall find 
below. 

Taking x\ — 0, x\ — r and x\ — t'>t> 0 the amplitude can be expressed as 

a . exp (~ik j± • x j± ) 

_ ^ k] x 

x exp [— ik z r' + ik\r' + ik z 2 (j' — r)] A (3.11) 


After some algebra we find, with an obvious notation for color traces, 



A = 

4m (q 

°) 3 /?t,l(A) |[Tr(caia 2 a 3 ) - 

Tr(aica 3 a 2 )] - 

i) 


+ 

[Tr(aica 2 a 3 ) - Tr(caia 3 a 2 )] 

ri)} 

(3.12) 

where 






Ax = 

l(Pi 

-q) 2 - 

- m 2 + ie][(pi -k 2 - k 3 ) 2 - 

- m 2 + ie][(pi - k 3 y 

! — m 2 + ie ] 

A 5 = 

l(Pi 

- h - 

k 2 ) 2 — m 2 + ie ] [(pi - k 2 ) 2 

- m 2 + ie][(p 1 - k 3/ 

) 2 — m 2 + ie] 

A 3 = 

l(Pi 

-h) 2 

-m 2 + ie][(pi -k 2 - k 3 ) 2 

— m 2 + ie ] [(pi -k 3 \ 

) 2 — m 2 + ie] 

A 7 = 

[(Pi 

-k 2 ) 2 

-m 2 + ie][(p 1 -h- k 3 ) 2 

- m 2 + ie][(p 1 - k 3/ 

) 2 — m 2 + ie] 


(3.13) 


R 0 g 7 ( d 3 ki d 3 k 2 -A- 
y/2Awm J (27t) 3 (27t) 3 



Pi i 

* n= 

<: - ^ P 


J/V 


"Pi 


^ k 1 A*2 A k 2 

a 1 , <3 2' a 3 i 

X * * 

1 2 3 


Figure 5: J/ip production amplitude induced by three scatterings. Eight diagrams 
are generated by attaching each Coulomb gluon to one of the quark lines. 
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In the g° —> oo limit we have 

(Pi ~ q) 2 ~ m2 — —2m 2 ± 

(Pi - kf ) 2 ~ m 2 + ie ~ q°[k- + k i± ■ (k j± + k k± )/q° + ie} 

(pi - ki - kj ) 2 - m 2 + ie ~ q°[k z + k z + k k ± ■ (ki± + k j± )/q° + ie](3.14) 

where i,j, k is any cyclic combination of 1, 2, 3. Integrating over k 2 and then over 
picks the poles 

k z 2 = k z - k{ + (£ 1± + fc 2± ) • £ 3± /g° + (3-15) 


and 


k{ — k“ + fcu_ • (/ c 2 jl + k%x)/Q° + 

= k z + fcu_ • (fc 2 _L + h±)/q° + 2k 2 ± • k 3 ±/q° + ie (3.16) 


As in the cases previously studied, the first longitudinal transfer k{ ~ k z puts the 
ccpair on-shell, while the succeeding transfers kf -C k z maintain this on-shcllness. 
This is because 


| ki± ■ kj± | -C m 2 x 


(3.17) 


holdsQ in both relevant ranges (|1.3| ) and (|1.4|) . We finally obtain 




x 

x 


vmq 


Rod 


7 


jnW 24 


d 2 ki d 2 k 2 


7T m 
3 




n 


rp a j 

1 A'. A j 
0 3 


exp(—ikj± ■ Xj±) 




J (2tt) 2 (2tt, 

/il_L * “1“ f'ca\d^a2CLzd ^2_L * ^3_L 


(3.18) 


From Eq. (|3.18|) we deduce the following: 


• The last momentum transfer fc 3 is hard. 

• One of the Erst two transfers is hard and the other is soft. 


Let us consider the two parts of the amplitude ( [3.181 ) separately, 

M = Mi + M 2 (3.19) 

2 Note that |fcij_ • kjj_\ <C mf also holds when one transverse momentum is of order 
m and the two others are smaller: \kj±\ <C m for j i. We have checked that this is indeed 
true in the calculation of {pfi) for transverse and also longitudinal J/k’ in the logarithmic 
approximation log(m 2 /Aq C£) ) 1 . 
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where the index indicates which is the first hard transfer. It is straightforward 
to see that the Ai 2 part factorizes in a way analogous^ to Eq. (|2.11|) , 


iM 2 


r j2u i 

= H J X)g(k 1 ) -> g(q + h, A')] ^ 

x iM[g(q + h, X')g(k 2 )g(k 3 ) -> J/il>(q + k)} 


(3.20) 


or in pictorial form 


soft hard hard 

* * * 

1 2 3 



(3.21) 

The soft transfer ki± thus mimics rescattering of the initial gluon. (As remarked 
above, however, the longitudinal part k\ is relatively big in the sense that it puts 
the incoming heavy quark pair on its mass-shell, as seen from Eq. (|3.16|) .) 


Factorization is somewhat more involved for the Adi term of Eq. (|3.19D , where 
Aqj_ is hard. A relation like 


J/y 


r <- 0 =£ 


d 2 b 


2±i 


_L( 1 ^<- 0 1 (I » - 


X' 


J/y' 


soft hard 

* * * 

1 2 3 


A' 


(2 ny 


(3.22) 

cannot hold, since the soft k 2 ± exchange sees a color which has been ‘rotated’ 
away from the incoming color c by the k\ gluon. Rather, the color structure of 
the XAi term in Eq. (|3.18|) can be pictorially represented as 


(iMi) color d C aidifcL2a3d 




1PJ3> 


soft hard 

* * * 




<d 


J/y' 

D= 


(3.23) 


'color 


'color 


On the other hand, forgetting the color factors and considering only the Lorentz 
structure of Adi, we easily see from Eq. (|3.18|) the validity of the factorization 


3 The difference in the expressions of m 2 j_ in Eq. ( |3.18| ) and in the J/tjj production amplitude 
appearing in Eq. ( [3.20 ), respectively to 2 + (fcu_ + k 2 ± + fc3_i_) 2 /4 and m 2 + ( k 2 ± + fc3_i_) 2 /4, is 
negligible. See the comments for the case of r/ c in the end of section 2. 
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formula 


(iM 


1 ) Lorentz 


= <? 




T~ 

hard soft hid 

* * * 

1 2 3 


j/y 

D= 


Lorentz 


= E 

A' ' 




2_L 


(2tt) 


x x \ 

xxx xs x x x \ 

<? ■ < ?+Jc 2 ) 


1 

2 g° 


<?+fc 2 



> 

T 


* 



Lorentz 


« * 

1 3 


I 


(3.24) 


Lorentz 


The differing structure of Eqs. (|3.23|) and (|3.24|) prevents us from factorizing the 
Mi part of the amplitude in a straightforward way analogous to Eq. (|3.20|) for 
M 2 - The soft gluon ( k 2 ) between the hard transfers (ki,k 3 ) changes the color 
structure of the hard production amplitude and, as we shall see below, ensures 
that there is no interference in the cross section between M\ and M 2 . It is 
convenient to summarize our result for the full amplitude (|3.18|) in the form 


M(4g^J/ip) = Mx + M 2 

3 

= n 

3 = 1 L 


1 A'. An 


fa2CL^d^ , caid^ / ^l “ 1 “ ^fca 1 dda 2 a 3 dM 2 j (3.25) 


wheref] M % = (Mi)Lorentz- For instance, Eq. (|3.24j) is explicitly 

“ d 2 k 2 ± . 


iM\ = E 


2 iMei[g(q, \)g(k 2 ) 
x iM [ g(q + k 2 , X')g(k 1 )g(k 3 ) - 


- (27T) 


J !i ) (q + k)\ 


g(q + k 2 , A')] 2^0 


(3.26) 


where 


M e i[g(q, X)g(k 2 ) -> g(q + k 2 , A')] 

= 2 q°g 2 [-e(A) ■ e(A')] 
M\g(q + k 2j X / )g(k 1 )g(k 3 ) -*■ J/i>(q + k) 


, exp (~ik 2± ■ x 2X _) 


k 2 


= —imq 


o 


Rog 5 


m\\p2A 


irm 


z Pt,l( X') 


x 


d 2 ki 


exp(-ik 1± ■ xi± - ik 3± ■ x 3 ±) 


ki± ■ k 3± 
k 2 h 2 

"'ll. 3± 


(3.27) 


4 We leave for convenience the -^= color factor for the bound state in the Lorentz part. 
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4 Cross sections and ^-broadening 


In this section we derive the expressions for the production cross sections and (p 2 ± ) 
of the rj c and J/-0 bound states, based on the amplitudes found in the previous 
sections. 


The differential cross section for scattering on static centres is 


d 2 a n 

d 2 P ± 


1 

( 2 tt ) 2 (2 q 0 ) 2 





(4.1) 


where A4 n is the bound state production amplitude induced by n scatterings 
(n > 1) and P± = 2pi±. In addition to the implicit color/spin summation and 
averaging, the average sign in Eq. (|4.1| ) denotes an average over positions of the 
scattering centres. 

In the case of one scattering centre, and assuming a uniform distribution of this 
centre with a constant volume density p, our description is equivalent to the 
standard infinite momentum frame description with a gluon distribution oc 1 /x 
(see Appendix A). For n > 1 centres we assume in the following that the centres 
are distributed independently with the same constant density p. Thus the average 
of Eq. ( |4.1|) is defined to be 

r. n _ 

(\M n \ 2 ) = I Y[(pd 3 Xi ) |Af„| 2 (4.2) 

i —1 

where the overline indicates the remaining color/spin sum and average. We also 
define the volume, transverse area and length of the target by 


V = J d 3 Xi 

S = j d 2 Xi± 

L = J dx\ (4.3) 


Our approach can in principle also deal with a case of correlated centres, which 
might be of phenomenological relevance, as discussed in the final section. 

From Eq. (|4~T|) we obtain the (P 2 ) induced by n scatterings as 


(p2) _ I d 2 P± P 2 d 2 a n /d 2 P ± 
± U f d 2 P±d 2 a n /d 2 P± 


(4.4) 


4.1 Tj c production 

One scattering 
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We first give the cross section for the basic single gluon exchange process. Squa¬ 
ring the amplitude M.(gg —■> r/ c ) of Eq. (|2Af) and using e = (0,1, ±i, 0)/\/2 for 
the projectile gluon polarization states we obtain 


dai_ = 4 tt 

dP 2 9 pv M(P]_ + M 2 ) 2 P]_ ’ 


(4.5) 


where M = 2m and the factor pE arises from the trivial averaging over Ti. 
Integrating over Pj yields 


4 tt w R 0 a* , 

(Ti = —pE ° log 
9 M 5 


M 2 


^QCD, 


(4.6) 


Two scatterings 

For the production process induced by two scatterings we use the convolution 
formula of Eq. (|2.11|) , which as we now shall show implies a similar convolution 
at the level of the cross section, 


d 2 a 2 

d 2 P ± 


[ d 2 k-\ i 


/ 

_d 2 fci_L_ 


J i 


d 2 cri(P± - fci±) 
d 2 P, 


(4.7) 


The first factor is the probability density for gluon elastic scattering off the first 
static centre. Thus d 2 cr 2 /d 2 P± is given by a convolution of the probability for 
having a gluon elastic scattering of transfer k±± first and the differential cross 
section for the basic production process with a single scattering of momentum 
transfer P± — k\± occurring thereafter. 

To derive Eq. ( |4.7|) we express the elastic gluon and p c production amplitude in 
the form 


iM e i[g(q, X)g{k 1 ) -> g(q + k u A')] = fca lc 'T^ iAl exp(-ik 1± ■ x 1± )e{\) • e(A') 

x iM et [q, fci] (4.8) 

iM[g(q + fci, X’)g{k 2 ) -»■ g c (q + k)] = T^ M exp (~ik 2 ± ■ X 2 x)e’ i (X') 

x iMfiq + h,P] (4.9) 


where the expressions for the A4 amplitudes are directly obtained from Eqs. (|2.8|) 
and ( CT - 

Squaring the amplitude ( |2.11| ) leads to 


d 2 <7 2 


d 2 k i i d 2 k' 


d 2 P 


± 


( 2 vr ) 2 J ( 2 tt ) 2 ( 2 ?r) 


i_L a 

-7 exp 


-i(ki± - k' 1L ) ■ (fi ± - x 2± ) 


x f m cf^T‘i A T< A . [E(A) ■ £ (A')] MX) ■ E(V)]- ( 4%M j (Ahthhl 


2 g c 


2 g° 
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X T%, M T% 2 A'^ (A ')e* (A") 


n'/\r\~»"r\"\* I + ki, P ] \ ( P] 


2 q° 


2 q° 


(4.10) 


Summing and averaging over colors, over positions x j, x 2 and over the interme¬ 
diate polarizations A', A" gives 


d2 ^ _ 1 />Z , 

1 pL 

Mei[q, ki\ 

2‘ 

1 pV 

e^M^[q + h,P] 2 

d?P L 2 J 

i 

to 

to 

2q° 


(2tt) 2 2N 

l 

O 

CN 


(4.11) 

The second factor in the integrand is readily checked to give the second factor in 


Eq. ([4.71) . The first factor is related to the probability density for gluon elastic 
scattering through 


d 2 Vel 

d 2 ki± 


1 d 2 a. 


el 


l l 


Sd 2 k 1± (27t) 2 5 


-A 4 el 


1 pL 


M 


el 


2 q° 


2 q° 

1 PL g 4 
( 2 tt ) 2 2 ki 


(4.12) 


i_l 


The overall factor | in Eq. ( |4.1ip arises from the constraint x\ < x 2 . The bracket 
notation [ ] i in Eq. ([4.71 ) is to be understood as implying that the gluon elastic 
scattering occurs before the hard transfer P± — ki±. Absorbing thus the factor | 
in the definition of the probability density, Eq. (14.71 ) follows from Eq. (|4.11|) . 


n scatterings 

The generalization of Eq. Q4.71) to the case of n scatterings may be written as 


d 2 


CT ri 


, 71—1 


d 2 P\ 


= n \d 2 k 


i± 


i— 1 


d 2 V, 


el 


d 2 k 


i± 


d 2 a 1 (P ± - E?=i kj 


■j-U 


d 2 Pi 


(4.13) 


Using Eqs. (ED, ED and (14.12|) we obtain the mean transverse momentum 
squared in n scatterings, 


Xl)n = 


. 71—1 

n 

i =1 


d 2 k, 


i± 


d 2 k. 


71 _L 


kt 


71 


iX. 


k\Mi . =1 


' 3 -L 


, 71—1 

n 

i=l 


d 2 kj 


i_L 


d 2 k 


(4.14) 


71 _L 


k il 




where the correct value M 2 to use is (see the discussion following Eq. Q2.26Q ) 

M\ = M 2 + Pl = M 2 + (k 1± + ... + k n± ) 2 (4.15) 
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To logarithmic accuracy, log(M 2 /Ag CD ) 1, the lower and upper limits of the 
logarithmic integrals in Eq. (|4.14|) can be set to A ‘q CD and M 2 , respectively. The 
result is 


(Pl)n = 


M 2 


log 


M 2 

A 2 QCD, 


+ in- l)K 2 QCD \og ) 

\ 1v qcdJ 


(4.16) 


We know from section 2 that in r] c production induced by n scatterings the last 
transfer k n _l is hard, whereas the (n — 1) first ones are soft. Thus we may write 
Eq. (|4.161) in a transparent way, denoting n so f t = n — 1, 


(Pl)n = £(*&} = «x> + n, c , t (kl, tL ) (4.17) 

2=1 


It is interesting to relate the mean value of n so f t to P e /. Integrating Eq. ( |4. 13| ) 
over Pj_ we get 

= [Pel] 1 • • • \Pel\n-X *1 = ^eT^l ( 4 - 18 ) 

The measured (P±) is obtained by summing over the number of rescatteriugs n, 


(Pi) 


En=l (PDnCTn ^n=l ^ 

E OO Y^°° rr 

n =1 a n 2 ^ n =1 


(P±)n =1 + (nsoft) {k SO ft ±) 


(4.19) 


For V e i C 1 we have 


(nsoft) = (n - 1) 


E oo 
n=2 


n 


1 )(T„ 


E oo 

n =1 a n 


*2 

01 



(4.20) 


Since (■ n so f t ) is given by the fraction of events where a single soft scattering 
precedes the hard one it is proportional to the length of the target. This soft 
scattering is responsible for the medium-induced P]_ of the r/ c , 

A Pi = \v„(kl, tl ) (4.21) 


which, according to Eq. (|4.16| ), depends on the mass M of the bound statef]. The 
derivation of the Pj_-broadcning of the J/i)j proceeds along the same lines, which 
we now summarize. 


5 We have assumed a fixed coupling constant a s 


just replaces the logarithmic dependence of the medium-induced ± j_ 
logarithmic dependence. 


Taking into account the running of a s 
Pf in Eq. (4.16) by a double 
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4.2 J/'il) production 


As we saw in section 3, when gluon radiation is neglected the basic J/A production 
process involves two hard scatterings. We quote only the transverse J/-0 cross 
section obtained from Eqs. (|3.10|) and (|4.1|), 


cr 2 = 


407T 

81 


M 7 


p 2 VL- 


log 2 


M 2 


A 2 

d v QCD, 


(4.22) 


The cross section for three scatterings is obtained from the amplitude given by 
Eq. (|3.25|) . Since d^dfca^ = 0, the interference between M. i, where ki± is hard, 
and At 2 , where k 2 ± is hard, vanishes and one gets a similar convolution as Eq. 


(14770 for the rj c , 


dV 3 

d 2 Pi 


1 

3 

1 

3 


d 2 k 


2_L 


d~k 


i_l 


d 2 V e i d 2 a\\P L - k 2± ) 

d 2 k 2 ± d 2 P± 

d 2 Vei d 2 qf(P L - ki±) 

d 2 P , d 2 P ± 


(4.23) 


The upper indices of a 2 indicate which transfers are hard and the factor 1/3 is 
due to the fact that the integration over longitudinal positions yields 1 / 2 ! for cr 2 
and 1/3! for < 73 . Absorbing this factor by defining the probability density to have 
the (soft) elastic scattering on a given centre i, 1 < i < 3, we find 


d 2 o-i 

d 2 P± 


= / d 2 k 



d 2 (? 2 3 (P± ~ k 2 ±) 


2 d 2 P ± 

d 2 af(P 1 - 


h±) 


d 2 k 1± \ 1 d 2 P ± 

Eq. ([4.24Q generalizes in a straightforward way to n scatterings, 


(4.24) 


d 2 a n 

d 2 P± 


n —1 r n —1 

E/n 


d 2 p 


j -L 


i= 1 


j'=l 


d 2 V, 


el 


d 2 k 


i -l 


d 2 a^(P x ^E]Zlk j± + k i± ) 


d 2 P 1 


(4.25) 


where i and n denote the hard transfers, while the n — 2 other transfers labelled 
by j are soft. Eq. (|4.25|) is valid for both transverse and longitudinal J/i]). 

The expression for (P±) n is then 

2 


(Pl)n = 


n 

i=l 


d 2 k; 


i L 


k A 


(fcn ■ k n ±)~ + ... + (fc n _ 



n 

i —1 


d 2 k 


z_L 




(k\± ■ k n ,j_) 2 + ... + (k n—l± ' k n ±) 




T,L 


J Ml 


(4.26) 
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where 


$ = 1 

AM 2 Pi _ AM 2 (^ r 
~ M^ _ n ' rA ' ^ ' 


Mf v • 1 

- 1 V=1 


i-L 


(4.27) 


In Eq. (|4.26|) it is necessary (for longitudinally polarized J/^’s) to keep the full 
expression for Ml which appears in Eq. ( 3.18|) . 

For transversally polarized J/^’s we hnd, similarly to the case of the r] c , 


< Pl)n = 2(kl ardX ) + (n-2)(k 2 soft± ) 

2 M 2 

3 


log 


M 2 

Aqcd, 


+ {n ~ 2) A 2 qcd log 


' M 2 ^ 

,^-Qcn, 


(4.28) 


Averaging over the number of scatterings n, and denoting n so f t = n — 2, 

(Pi) = (Pl)n=2 + (n soft )(k 2 soft± ) (4.29) 

From Eq. (^4.25|) we have 


On. — 2 . P e I <72 


n\ 


(4.30) 


which for V e i -C 1 yields 


(n S oft) — (n 2) 


E”. 3 (n ^ 2K 

E oo 

n =2 a n 


O 3 

02 


V e , < 1 


(4.31) 


The medium-induced (Pi) for transversally polarized J/-0 is thus 

API = \v a (kl, tl ) (4.32) 

and depends logarithmically on M. We note that | V e i is the probability to have 
one soft gluon scattering before the second (and last) hard transfer. 

In the case of longitudinally polarized J/ij) the expression for (Pi) is more in¬ 
volved since the second term in 


<Pl>» = £<&> + 2 £<*■-i-M (4.33) 

i i<j 

is non-vanishing. However, we hnd that API still depends logarithmically on M. 
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5 Discussion 


In the present work we studied rescattering effects in hard collisions. Most PQCD 
calculations so far have concentrated on the hard vertex itself, thus involving 
only a single projectile and target parton. We were particularly motivated by the 
discrepancies observed between theory and data for quarkonium hadroproduction, 
and the evidence for large nuclear effects both in the cross section and in the 
average transverse momentum. 

To facilitate the calculations we assumed the target scattering to occur off fixed 
centres, which selects Coulomb exchange. As illustrated in the Appendix, at least 
for a single scattering this is equivalent, in the high energy limit ( x tar g —> 0 ), to 
a standard calculation with a specific target gluon structure function. Since the 
assumption of Coulomb exchange simplifies the calculations considerably, it will 
be worthwhile to investigate how general this equivalence in fact is. 

As we already discussed in the Introduction, our calculation shows how soft scat¬ 
tering factorizes from the hard vertex. This is particularly delicate in the case of 
J/ip production, where soft scattering may occur between the two hard exchanges, 
and thus affect the color structure of the hard vertex itself. An unexpected result 
of the calculation was that the (relatively) large longitudinal momentum transfer 
required to put the heavy quark pair on its mass shell always occurs at the first 
scattering, even if that scattering is soft in terms of the transverse momentum 
exchanged. This implies that our analysis is applicable to fixed target data on 
T production, even though the beam energy is not large enough for the bb pair 
to travel very far off its mass-shell. Note that these general features remain true 
when considering J/ip or T photoproduction in our model. 

The exchange of two hard gluons, required for radiationless J/if) production, is 
typically a higher twist effect and thus is suppressed by a power of the heavy 
quark mass mg. This is because of the small probability to find two partons 
within a transverse distance of order 1/mg from the heavy quarks. However, our 
calculation focusses attention on the fact that there actually are partons within 
such close distance, namely those created by the Q 2 evolution of the projectile 
and target partons. Thus the projectile gluon of Fig. [1] in reality is not on-shell, 
but has a logarithmically distributed virtuality due to previously emitted gluons. 
It is ‘hard’ in the sense of Eq. (0), and due to its relatively short life-time must 
still be at a short transverse distance from its radiated partners. Could gluon 
exchange between the quark pair and such partners be an important effect in 
J/'ijj formation? 

Such a contribution is in fact nothing but a higher order loop effect in the J/^j 
production amplitude, and hence is definitely of leading twist. It is suppressed by 
a power of a s in the cross section compared to the lowest order process involving 
gluon emission. However, this is at least partly compensated by the fact that 
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no energy is lost from the quark pair. Such exchanges could thus be particularly 
important for the Tevatron data on J/ip production at large p±, due to the strong 
‘trigger bias’ effect which favors hard fragmentation mechanisms. Although a full 
higher order calculation of J/ip hadroproduction is a challenge that remains to be 
met, the particular contribution from rescattering off evolutionary gluons should 
not be difficult to estimate. 
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Appendix 


A Relating the heavy target quark density to 
the gluon distribution 


In this Appendix we wish to demonstrate the equivalence between our heavy 
quark target model and the standard structure function formulation. In parti¬ 
cular, we derive the relation between the number of heavy quarks and the gluon 
distribution G(x) in the limit x —> 0. We choose to consider the process yg —> qq 
(see Fig. |6|), where massless quarks are produced at large p±. We expect the 
equivalence to hold similarly for any other hard process in the small x limit. 



Pi 

P 2 

t 



Figure 6: Amplitude for the process yg —*• qq. 

The standard infinite momentum frame expression for the y g —> qq cross section 
in the limit where the photon energy v —> oo while the transverse momentum p± 
of the produced quarks and their squared invariant mass s — xs — 2xumN are 
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kept fixed is 


da 

dp\ dz 

where p\ = z( 1 — z)s. 


ne 2 aa s 

Q b 

p\ 


xG(x) z 2 + (1 



(A.l) 


In the target rest frame (see Fig. [|) z is the fraction of the photon longitudinal 
momentum carried by the quark, 


Pi = (Pi±, zv ) 

P2 = (P2±, (1 - z)v + k z ) (A.2) 

2 vk z = -PiJz-p 2 2 J(l - z) 


The amplitude for photon scattering on a heavy quark t in the limit we consider 
can then be written 0 


Ti. 'tt -» qqt) = ee,g Ta B T£ d 


Amtvz(l — z) 


X 


u(pi)e- 7 v(jp 2 - k) u{pi - k)e- ^fv(p 2 ) 


pi± 


p'L 


(A.3) 


After summing and averaging over spins the square of the amplitude simplifies 
to 

,N 2 -lAmWzil- z) r 

I Z4-7T I * 


\T\ 2 = 16e^a^(47r) 3 ' 


8 A" k 2 ± p 2 ± p 2 x 
In the expression for the cross section 


z 2 + {l-zf 


(A.4) 


da 


d 2 k i 


dp\ L dz 4m^(47r) 4 J z(l — z)v 2 


\T\‘ 


(A.5) 


the /cx-dependence of the integrand is (for k± -C pu, p 2 ±) of the form 1 /k\. 
Here we introduce R ~ Rqcd as an infrared cut-off for the logarithmic integral, 
where R is the nucleon radius, playing the role of a color screening length. To 
leading logarithmic accuracy the scattering occurs incoherently over the pointlike 
target quarks. The upper limit of the /ex-integration is set by k± ~ pix, where the 
approximation p^p = (k± — pi ±) 2 — p 2 ±_ fails and the integrand is more strongly 
damped in k±. This defines the logarithmic approximation as 


log 


P± 

A qcd 


> 1 


(A.6) 


In order to deal with measurable cross sections one has to average over the possible 
positions of the heavy target quark. We assume in this paper a uniform density 
for the target quark, denoted by p, and thus we multiply Eq. (|A.5|) by a factor 


26 























pV counting the number of static centres in one hard scattering, V being the 
target volume. 


Comparing the resulting expression for the h^-integrated cross section with that 
of the infinite momentum frame result (|A.1|) we find 


xG{x)=pV ^ir aM)log {i^) • (A - 7) 

or, taking N = 3, 

pV = 33 ~J nf xG(x ) (A.8) 

16 

We have checked that the relation (|A.8|) is obtained also if one considers r) c 
production through two-gluon fusion, by simply comparing the expressions for 
the cross section in the target rest frame obtained from Eq. (|2.8|) with that of the 
infinite momentum frame (see, eg, Ref. [[ZSj). 
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